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Context. Z CMa is a binary system which consists of two young stars: A Herbig AeBe component "ZCMaNW" embedded in a dust 
cocoon and a less massive component "ZCMaSE", which is classified as a FU Orionis type star. Associated to the binary system 
is a giant parsec size jet. Past spectropolarimetric observations showed that the position angle of the linear optical polarization is 
(~| I perpendicular to the jet axis, indicating that the visual light escapes the cocoon via cavities aligned with the jet axis and is then 

scattered back into the line of sight of the observer Recently, the system showed the largest outburst reported during the almost 90 
I ■ years of available observations. 

O ' Aims. We present new spectrophotometric and spectropolarimetric data obtained in 2008 during the recent outburst phase. 

Methods. The data obtained in the visual spectral range at medium spectral resolution were used to study the geometry of the system 
C/3 , from the linear polarization spectra and its magnetic field from the circular polarization spectra. 

Results. During the recent outburst we detect that the ZCMa system is polarized by 2.6% in the continuum and emission line 
spectrum, with a position angle still perpendicular to the jet. From the high level of polarization we conclude that the outburst is 
associated with the dust embedded Herbig AeBe NW component. The deep absorption components of the Balmer lines in the velocity 
frame extending from zero velocity reaching a wind velocity of ~700 km s"' , together with the absence of a red-shifted broad emission 
at similar velocities, indicate a bi-polar wind. We do not detect a significant mean longitudinal magnetic field during the outburst, but 
• the data obtained in 2004 indicates the possible presence of a rather strong magnetic field of the order of ~1 kG. However, we 
' critically review the applied method of magnetic field measurements in the presence of a strong stellar wind. The main result 
' of our studies is that the bolometric luminosity of Z CMa remained surprisingly constant during the recent "outburst". We 
conclude that either the geometry of the cavity through which the light escapes from the cocoon has opened a new path, or 
' that the screen of dust, which reflects the light toward the observer became more efficient causing the observed increase of the 
visual brightness by about 2? 5. 

' Key words. Stars: individual: ZCMa — Stars: pre-main sequence, winds, outflows, magnetic fields, variable: general, binaries: close 

1. Introduction (119951) resolved the binary system with speckle interferometry 

' at 730 nm and 656 nm. From the data obtained in the quiescent 

■ ZCMa (HD 51370) is a young visual binary system consisting ^^^^^ 1989.84 the authors deduced L = 1300 Lo, R ^ UR^ 

■ - - ' of a 1 2 Mo Herbig AeBe star and a 3 Mq FU Ori companion at ^n^j j^^.^ ^ qOOO K for the Z CMa SE FU Ori component and L = 

an age of 0.3 Myr (e.g. Alonso-Albi et al. |2009D. Both com- 2400 Lq, R = 1690 /^o, and Tbb = 980K for the near infrared 

ponents are surrounded by active accreting disks. Associated zCMaNW component. The authors attiibuted the excess flux 

wi th the bmary system is a giant parsec size jet (Poetzel et associated with ZCMaNW at optical wavelengths to scattered 

al. II989II. First evidence for the binanty was an elongation of jig^t, arguing that photons originating from the hidden central 

the K-band image of this star toward the northwest (Leinert & go^j-ce are scattered back into the line of sight. Scattering of the 

Haas 1987 and Koresko et al. 1989 ). The binary was resolved photons is assumed to take place along the walls of a cavity (usu- 

by Koresko et al. (1991) with near-IR speckle interferometry. ally ascribed to a sti'ong bipolar outflow) evacuated within the 

The angular separation of the binary was only -^.'l with a po- envelope. Radio observations in 2005 and 2006 in the cm and 

sition angle of 120°. At that time the south-eastern (ZCMaSE) ^m range suggest that we see ZCMaNW at an inclination of 

component was dominant in the visual spectral range and near- 30° (Alonso-Albi et al. 2009) . This means that we do not observe 

infrared J and H bands, while the north-western component was t^e inner walls of the bipolar cavity within the dusty envelope, 

found brighter at 2.2 yum and longer wavelengths. Thiebaut et al. vvhere the photons are scattered 



Send offprint requests to: Thomas Szeifert The scattered light scenario is weU consistent with optical 

* Based on observations made with ESO Telescopes at the Paranal spectropolarimetric observations of Whitney et al. ( 1993), who 

Observatories under programme IDs: 074.C-0442, 081.C-0410, and detected in 1991-1992 a significant linear polarization of 1% 

282.C-5041. to 2% in the continuum and up to 6% in the strong emission 



2 



Szeifert et al.: The nature of the recent extreme outburst of the Herbig Be/FU Ori binary Z CMa 




11 - 

11.5 E \ < \ < Li 

1990 2000 2010 



Fig. 1. Light curve of ZCMa over the last decades. Open sym- 
bols denote the observations from American Associations of 
Variable Stars Observers (AAVSO). The filled symbols corre- 
spond largely to the ASAS data and during the peak of the 
lightcurve to our measurements (blue asterisks). Epochs of 
NACO observations, the high resolution images as summarized 
by Thiebaut et al. ( 119951 ) and the epoch of spectropolarimetric 
observations (arrows) are indicated on the top. 



lines. The different polarization of continuum and lines have led 
Whitney et al. ( |1993l l to suggest that the emission lines have their 
origin in the Z CMa NW component. At that time the NW com- 
ponent contributed about 20% of the flux of the whole system 
at visible wavelength. The Z CMa SE component is assumed to 
be unpolarized, the NW component is assumed to show 6% lin- 
ear polarization and the composite spectrum of the binary would 
then depend on the flux ratio at the respective wavelength. The 
position angle of the polarization of 150° in November 1991 and 
154° in March 1992 was roughly perpendicular to the 3.6 par- 
sec scale jet of ZCMa at a position angle of 60°, discussed by 
Poetzel et al. ( I1989I I. The same authors reported several velocity 
components in the jet of which the fastest knots reaches a ve- 
locity of about -620km s"'. This position angle of polarization 
indicates that the light from the Z CMa NW component primar- 
ily escapes from the envelope by scattering off the walls of a 
jet-blown cavity. A smaller scale jet with a similar position an- 
gle was recently discovered in X-rays by Stelzer et al. (2009). 

According to photometric and visual observations available 
from various archives, Z CMa exhibited recently the largest out- 
burst reported during the almost 90 years of available observa- 
tions. The light-curves presented in Fig. [1] were extracted from 
the archives of the 'All Sky Automated Survey" (ASAS) re- 
jecting the data points with quality flags C and D (Pojmanski 
I2002I I. More data over a larger time span were retrieved from the 
AAVSO data base. For the later visual estimates the data points 
taken for a given month were averaged after clipping data points 
beyond a 2cr threshold. Irregular variability was observed in the 
thirties of the last century. Later Z CMa steadily increased its 
brightness, reaching a phase of constant brightness until about 
1973. Since then the star was fading from almost 9"' in 1973 
to 1CP5 in 2007. Superimposed to the fading lightcurve we can 
identify in Fig. [1] several more recent outbursts, the most promi- 
nent ones were observed in 1987 and 2000 with a duration sim- 
ilar to that of eruptions of EX Lupi-type (EXors) young stars. 
Between 2002 and 2006 the light-curve appeared very irregular 
A gradual increase of the visual brightness started at the end of 
2007, reaching the highest brightness between November 2008 
and May 2009. This bright phase has now ended, with the latest 
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Fig. 2. Spectrophotometric flux of Z CMa in 2004 (lower spec- 
trum) and 2008 (upper spectrum). 
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Fig. 3. Spectropolarimetry of Z CMa in Hy. 

observations reporting a visual brightness fainter than 10™ after 
a very rapid decline. 

In this letter we present new optical spectropolarimetric ob- 
servations of linear and circular polarization in the Z CMa sys- 
tem, obtained during the recent outburst, and discuss the cause 
for this outburst. 

2. Results and Discussion 

2. 1 . Spectroscopy and Linear Polarization 

The spectropolarimetric data sets provide very important in- 
sights into the nature of the at visual wavelengths unresolved 
binary star We found significant changes of the line shapes in 
the linearly polarized light with respect to the data published 
by Whitney et al. ( I1993I ). In our data obtained during the ex- 
treme outburst, we see a flat polarization spectrum with an av- 
erage polarization P=2.6% +0.1% at 4350A at a position angle 
of 160° + 1°. This polarization angle is roughly perpendicular 
to the parsec scale jet associated with Z CMa, and similar to ob- 
servations presented by Whitney et al. ( 119931 ) who measured in 
November 1991 P = 2.0%, PA = 150° + 1.5° and four months 
later P - 1.4%, PA = 154° + 1.5°. However, in the deep absorp- 
tion troughs of the Balmer lines and other strong wind absorp- 
tion lines, we detect depolarization signatures (see Fig.|3]l. These 
line components caused by the accelerated wind show a com- 
posite spectrum: We observe a broad blue trough typically as- 
sociated with an accelerated wind, visible to about -600 km s"', 
with less deep absorption up to about -800 km s"^ As is shown 
in more detail in Fig.|4] this is not mirrored by emission on the 
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Fig. 4. Logarithmic spectral shape of the Balmer lines Hq, to H9. 
The variability is most pronounced in the absorption troughs and 
forHo, in the symmetric electron scattered broad wings. 



red side of the lines as would be expected for spherical sym- 
metric winds. The trough is rather deep, indicating that a large 
fraction of the stellar disk is obscured by the wind accelera- 
tion zone. Superimposed to the wind feature is a very narrow 
emission component of 50 to 60kms ' half width, which can 
be seen as well for faint low excitation lines like of Fe 11 and 
similar species. In contrast to the observations of Whitney et al. 
( I1993I I. we find that the narrow emission lines show the same de- 
gree (2.6%) and position angle of polarization as the continuum. 
These observations can be understood within the above men- 
tioned model that the visible light can only escape via cavities 
from the dust cocoon, most likely aligned with the bi-polar jet 
and a changing flux ratio between the two stars: In our recent ob- 
servations the continuum radiation as well as the emission com- 
ponents are dominated by the NW component, while only in the 
deep absorption troughs the unpolarized SE component is still 
contributing. Whitney et al. ( 1993) observed the binary in a dif- 
ferent state, when the unpolarized SE component dominated the 
flux in the near-UV and visual. The strongest emission lines pri- 
marily emitted by the NW-component reached up to a polariza- 
tion of about 6%. This high value must be seen as the degree of 
polarization of the NW-component at all wavelength and accord- 
ingly the polarization of the NW-component was higher at that 
time than what we see now in the visual and near-UV continuum. 
The change of the degree of polarization of the NW component 
from 6% at a low state to 2.6% at the brightest state, together 
with the brightening of the target indicates that the geometry of 
the cavity through which the light escapes from the cocoon has 
opened a new path, or that the screen of dust, which reflects the 
light toward the observer, became more efficient, causing an in- 
crease in the visual brightness by about 2™5. In the strongest 
Balmer lines the absorption trough presents a double plateau 
with the deepest depolarization starting at the terminal wind ve- 
locity of about -600 km s ' and ending at about -200 kms '. 
At red-shifted wavelengths the NW component dominates the 
composite spectrum. It is possible to explain the double plateau 
on the blue-shifted side, if the spectrum of the SE component is 
that of a fast rotator with hydrogen absorption lines such that the 
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Fig. 5. Circular polarization in the Balmer line spectra observed 
in November 2004 (on the bottom) together with the spectral 
gradients dl/dA. The spectra are plotted in the velocity frame to 
display Balmer lines from 4861 A to 3734.370 A. Both curves 
indicate that the s-shaped profiles of the high-n Balmer lines are 
not centered at the rest frame wavelengths. For the low-n Balmer 
lines, in this figure strongly deviating in dl/dA from the high-n 
profiles, the peak of the emission is close to the rest frame wave- 
length. Accordingly the spectral gradient is zero at the rest frame 
wavelength and doesn't correlate with the Stokes V profile. 



contribution of the SE component to the composite spectrum is 
smaller from -200 to km s"' at wavelengths. Other line species 
which show depolarization signals in the blue-shifted parts of 
line profiles are calcium lines in which the depolarization signal 
can only be seen at v < -200 km s ' and the Fe ii-lines, which do 
not present the deep absorption troughs in the flux spectrum, but 
show a double absorption feature in the linear polarization with 
an absorption at the velocity of the fast wind and another absorp- 
tion at the system velocity. For the later we can not conclude if 
the depolarization has an intrinsic origin in the NW component, 
or if the presumed young FUOri star contributes unpolarized 
line emission to the composite spectrum. 

We note that the depth and the width of the absorption at 
the blue-shifted side of the Balmer lines, in absence of a broad 
600 km s"' red-shifted emission, indicates a bipolar geometry of 
the wind. If no such broad red emission is observed, then it can 
be argued that either there is no gas streaming out except out 
of the polar region, or that the recombination radiation in the 
line is not visible through the cavities, but blocked by the disk 
or cocoon. For the Ha line we find however a broad emission 
component beyond the velocity of the wind absorption feature. 
The Ha absorption trough is not as deep as for other Balmer 
lines. This broad component is caused by electron scattering as 
is often seen in luminous hot stars and therefore not related to 
high velocity outflows. 

The spectral slope as indicated in Fig. |2] appears redder dur- 
ing the recent outburst, which contradicts an increase of temper- 
ature as source of the luminosity change as proposed by van den 
Ancker et al. (2004). Similarly, the very insignificant change of 
the depth of the Balmer jump would argue against a dramatic 
change of stellar parameters. 

2.2. Circular Polarization 

The analysis of circular polarization measurements obtained 
during the outburst reveals no significant detection. In the data 
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obtained in May 2008 we detect a weak signal of circular po- 
larization in the trough of all Balmer lines, which however does 
not match in wavelength the unpolarized line profile. However, 
while reviewing the archive data akeady published by Wade et 
al. (120071 we do detect a clear circular polarization signal in par- 
ticular at the blue shifted side of the high-n Balmer lines. For the 
low excitation lines from HyS to He we measure -160 + 84 G, 
while from the higher excitation Balmer lines a very significant 
field of -1231 + 164 G (see Figs. Hand HI was obtained. Still, 
these results have to be taken with some caution. To properly 
interpret the circular polarization signal, it is necessary to con- 
sider the composite nature of the binary system spectrum, the 
superimposed effects of the stellar wind, and the variability of 
the target. The B-field diagnostics as applied above correlates 
the s-shaped profile in VII with the s-shaped profile of the spec- 
tral gradient d//d/l of the Voigt profiles. The B-field diagnostics 
may, for instance, be disturbed by the light of the companions, 
and, in particular, by the stellar wind features due to the spectral 
gradient d//d/l along the x-axis of the diagnostic diagram (see 
Figs. |5]and |6]l. It should as well be considered that the target 
is variable on very short time scales, such that spectra may have 
changed between the respective exposures at +45°. 

It is as well unclear if these fields were detected in the dust 
embedded Herbig Ae star or in the SE companion. In the first 
case the strong field is only seen in the fainter Balmer lines, 
because the high opacity of the expanding wind at the blue 
shifted part of the strong lines would not allow light from the 
photosphere to escape from the NW component. Alternatively, 
without any indication of the contribution of the respective two 
stars to the UV and visible spectra, one can as well explain the 
data if the magnetic field was detected in the SE component. 
This component may have dominated the spectra of 2004 in the 
wavelength range of the Balmer lines where we measure the 
Stokes V features, whereas later in May 2008, the contribution 
of Z CMa SE became insignificant on the blue-shifted side of the 
P Cygni line of Z CMa NW such that we detect a hint of circular 
polarization in the red-shifted part of the hydrogen lines only, 
but no S-shaped feature, which would lead to a magnetic field 
detection. Nor would the wind dominated spectral line gradients 
d//d/l linearly correlate with Vjl. With the few available spec- 
tra it is not yet possible to reliably conclude on the origin of the 
features in the circular polarization spectrum. 

2.3. Origin of tine Outbursts 

Despite of the different origin of the NIR Ks band flux from 
the surrounding dust cocoon and the optical flux with an origin 
from the pole of the star, we presume that a strongly changed 
optical/UV source would heat up the dust and the increase of 
the luminosity of the stellar core would be re-emitted by thermal 
radiation in the Ks-band and beyond. There is however no sig- 
nificant change of the Ks-band flux as far as we can tell from the 
K-band flux ratios given in Table[2]or from comparing the most 
recent report of infrared observations by Antoniucci et al. ( I2009I I 
of Ks= 3.56 at the peak of the lightcurve in March 2009, with 
Ks= 3.79 reported by Koresko et al. (1991). Similarly there is no 
change in the X-ray properties reported by Stelzer et al. ( I2009I I 
in respect to earlier measurement at a 'quiescent' phase. If the 
short-term optical brightening was an EXor-like accretion event. 
X-ray variability would have been expected. Our polarimetric 
data leads to a new interpretation in which the X-ray emission 
associated with the binary may be attributed to the lower-mass 
SE component which is not involved in the outburst. 



The apparently insignificant change of the K-band thermal 
emission rather argues against the strongly increased luminos- 
ity of 3 X 10^ Lq during short outbursts, as estimated by van 
den Ancker et al. (2004) from the visual SEDs with respect to 
only L - 2400 Lq estimated by Thiebaut et al. (1 19951 ) from the 
SED including the thermal emission. It rather underlines our in- 
terpretation that the bolometric luminosity during the outbursts 
remains constant, but the amount of visible light transmitted 
through the envelope has changed. The derived high luminos- 
ity was primarily a consequence of the early spectral type es- 
timated for the NW source from the He i and O i lines, which 
lead to very high bolometric and extinction corrections. We de- 
tect these lines with blue shifted broad absorption profiles and 
with P Cygni profiles respectively, both indicating an origin of 
the lines in the accelerated wind rather than in the photosphere 
of the NW component. The wind however, with indication of 
bipolarity and eventually with indications for a magnetic field, 
may well be excited by other sources rather than by the thermal 
radiation from the star. The short-term outbursts are neverthe- 
less superimposed to long-term FU Ori like changes in the vi- 
sual brightness, with a slow brightening by one magnitude until 
about 1970 and a subsequence fading by a similar amount. It is 
well possible that any of the binary components contributes to 
these changes. 

3. Summary 

We have shown that the NW component of Z CMa was undergo- 
ing a historically bright outburst in 2008/2009. Our polarimetric 
data supports a variable scattering geometry of the NW enve- 
lope as the cause for the optical outburst, rather than FU Ori or 
EXor events from the SE component. From the changes in the 
linear polarization spectra taken at the outburst phase we con- 
clude that the scattering path in which the light of Z CMa NW 
escapes from the dust cocoon must have changed. The linear 
polarization from the now dominant NW component was found 
smaller than the 6% estimated at the fainter phase in 1991. We 
do not find changes in the K-band magnitude in contrast to the 
recent change of 2™5 in the visual. A significant fraction of the 
variability may be explained by the changed, more direct line of 
sight in which we now see the Z CMa NW: More light from the 
central source would reach the observer through a widened cav- 
ity, while the central source directly observed in the K-band did 
not significantly change in luminosity. 
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Table .1. Journal of spectroscopic observations. q.oz | — , — , — ^ — ^ 



Date 


lil^Li U-lilCllL 


A range [A] 


Tlf*tf»r'tr»i' frnctn 
l-'CLCL.LUi, gllaiii 


2004/11/22 


FORS 1 


3470 


-5870 


TEK, 600B 


2008/05/22 


FORS 1 


3310 


-6200 


E2V, 600B 


2008/12/26 


FORS 1 


3310 


-6200 


E2V, 600B 


2008/12/26 


FORS 1 


3680 


-5120 


E2V, 1200B 


2009/02/23 


FEROS 


3570 


-9210 


Echelle 


2009/03/13 


FEROS 


3570 


-9210 


Echelle 


2009/03/26 


FEROS 


3570 


-9210 


Echelle 



Bz = -1231 +- 164 

MJD = 53331.307 ; hi^ih n Balmer line 




Table .2. Journal of NACO observations. The position angle is 
given with respect to the SE component. For the L-band data the 
angular separation and position angle was set to a fixed value to 
obtain a reliable flux ratio. 



MJD = 53331.307 ; 
Bz = -160 +- 84 



Balmer line 



Year 


Band 


fjvw 


/fsE 


PAn 


d["] 


2002/12/05 


L 


49.57 


+0.15 






2002/12/05 


J 


0.94 


±0.01 


309.90 


0.107 


2002/12/05 


Ks 


3.97 


±0.01 


311.9 


0.102 


2005/02/16 


Ks 


4.7 : 


tO.Ol 


311.43 


0.098 


2005/02/16 


H 


2.79 


±0.01 


310.98 


0.114 


2005/02/16 


NB2.17 


4.34 


±0.05 


314.74 


0.104 


2006/03/13 


NB2.12 


3.63 


±0.01 


311.93 


0.101 


2006/03/13 


NB2.17 


3.42 


±0.01 


312.43 


0.101 



Appendix A: Observations and Data Reduction 

The observational material discussed in this section includes 
FORS 1 spectropolarimetric data (Appenzeller et al. 1998) ob- 
tained in the framework of our programs carried out in May 
2008 and in December 2008. Additionally, we used ESO archive 
data from 2004 obtained in circular spectropolarimetric mode. 
These data were used by Wade et al. (I2007I I to search for the 
presence of longitudinal magnetic fields in a larger number of 
Herbig AeBe stars. These older observations (see Table [Tl) were 
carried out with a slit width of 0'.'8 and grism 600B (R » 1000). 
The flux spectra in the ordinary and extraordinary ray of the 
Wollaston prism were extracted for all position angle settings 
of the retarder waveplate. Since there were no available flux cal- 
ibrators observed for this setup in the respective night, we used 
flux calibrators taken in normal spectroscopic mode without the 
Wollaston prism. The new spectropolarimetric data from 2008 
were obtained with the same instrument, but now equipped with 
a higher UV response detector system with a smaller pixel size. 
During the December run we have used the higher dispersion 
grism 1200B in addition to the data with the 600B grism. In 
all our observations the slit width was set to 0'.'4 to optimize 
the resolution for the polarimetric measurements (R x 4000 and 
R =s 2000, respectively). The spectrophotometric calibrators re- 
trieved from the archive were observed without polarization op- 
tics on different nights. One source of uncertainty in the flux 
calibration is therefore the unclear spectral response of the wave 
retarder plates, for which we presume a high throughput over the 
full spectral range. However, significant slit losses are expected 
with this setup due to the narrow slit widths and the difficult ac- 
quisition of the target on the slit considering the extreme bright- 
ness of the target for 8 m-class telescopes. To remedy this dif- 
ficulty related to the absolute calibration, we re-scaled the flux 
at 5500 A to the monthly average of the AAVSO visual obser- 
vations of 9^4 in November 2004 and 8.4 in December 2008. 
In Fig. |2]we present spectrophotometric fluxes for the ZCMa 
system in 2004 and 2008, respectively. Further, for the proper 
interpretation of the measured polarization angle, we corrected 



-4.0 10-6 -3 10-& 0.0 lOO 3.0 10-6 4,0 10-6 

-4.67 lO.Qe-13 ^**3/r dl/dA 

Fig. 6. Magnetic field measurements using the data obtained in 
November 2004. The longitudinal magnetic field is derived fol- 
lowing the equation V/I - -geffC,/l^ from the slope of 
the linear fit through the data points (see Bagnulo et al. 2002} . 
Large open and small filled symbols are obtained from the high 
number and low number Balmer lines, respectively. 



the result with the chromatic zero angle of the super-achromatic 
half wave retarder plateQ To perform linear polarization mea- 
surements, a Wollaston prism and a half-wave retarder plate ro- 
tated in 22.5° steps between and 67.5° were used. For circular 
polarization measurements the quarter-wave retarder plate was 
used at the positions -1-45° and -45°. Our FORS 1 data were sup- 
plemented by three high-resolution FEROS spectra {R-42, 000) 
obtained in February and March 2009 to study the behavior of 
hydrogen lines. The journal of spectroscopic observations is pre- 
sented in Table [Tl 

To better understand the parameters of the resolved binary 
system we retrieved all available NACO adaptive optics near- 
infrared data (Lenzen et al. 2003) from the ESO archive. As for 
the FORS 1 data sets, the target is rather bright for imaging mode 
observations with 8 m-class telescopes and accordingly the vari- 
ous data sets were taken with a mixed choice of neutral density 
and narrow band filters. As a consequence we cannot provide 
absolute photometry for NACO data with the required level of 
confidence, but only the flux ratio of the two components (see 
Table O and the astrometric parameters of the binary. In case of 
normal jittered observation sequences we have used the Eclipse 
data reduction package provided by ESO to NACO users. In 
other cases we had only two acquisition exposures or a few jitter 
positions. 



' The data is available on the instrument web pages: 
http://www.eso.org/sci/facilities/paranal/instruments/fors/,. 



